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We recall the Einstein velocity addition on the open unit ball B of M.^ and its algebraic 
structure, called the Einstein gyrogroup. We establish an isomorphism between the Einstein 
gyrogroup on B and the set of all qubit density matrices representing mixed states endowed 
with an appropriate addition. Our main result establishes a relation between the trace 
metric for the qubit density matrices and the rapidity metric on the Einstein gyrogroup B. 
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1. Introduction 

Einstein addition is the standard velocity addition of relativistically admissible vectors 
that Einstein introduced in his 1905 paper that founded the special theory of relativity. In 
his book [3J the presentation of Einstein's special theory of relativity is based on Einstein 
velocity addition law. It also allows the reader to take the properties of the underlying 
hyperbolic geometry. 

A. A. Ungar has introduced and studied in [5] gyrogroups that are algebraic settings of 
hyperbolic geometry, and suggested three examples of gyrogroups corresponding to three 
models of hyperbolic geometry. It has been known that gyrogroups correspond equivalently 
to loop structures; see [3]. In Section 2 we review the Einstein gyrogroup on the open unit 
ball of the Euclidean three-dimensional space corresponding to the Beltrami-Klein ball 
model of hyperbolic geometry. 

Bloch vectors in the unit ball of M'^ are well-known in quantum information and compu- 
tation theory. A qubit density matrix, a 2-by-2 positive semidefinite Hermitian matrix with 
trace 1, is completely described by the Bloch vectors. Chen and Ungar have computed the 
Bures fidelity of qubit density matrices generated by two Bloch vectors and showed some 
equivalent formulas in their papers [6] and [7]. The Bures fidelity that plays an important 
role for the geometry of quantum state space has been manipulated into a form possessing 
distance properties, called the Bures metric. On the other hand, we apply the trace metric 
to the space of all invertible qubit density matrices, investigate its properties, and show in 
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Section 5 the relation between the rapidity metric of the Einstein gyrogroup on B and the 
trace metric of the qubit density matrices. 

2. Einstein addition and gyrogroups 

The velocities in Einstein's theory of special relativity are three-dimensional vectors with 
magnitude bounded by the speed of light. We assume the speed of light is normalized by the 
value 1, and call such velocities admissible vectors. The relativistic sum of two admissible 
vectors u and v of norm less than 1 is given by 

u e V = — ju + + ^!^-(u^v)u| , (2.1) 
1 + u-* V [ 7u 1 + 7u J 

where 7u is the well-known Lorentz factor 

^" = 7™- 

Note that u-^v is just the Euclidean inner product of u and v written in matrix form. 



Definition 2.1. The formula (j2.ip defines a binary operation called Einstein addition on 
the open unit ball B = {u : ||u|| < 1} of M^. 

The Einstein addition can be naturally defined on the open unit ball B of n-dimensional 
space M". See [21 Chapter 1] for a derivation of the Einstein addition law. 

Remark 2.2. The Einstein vector addition on the open unit ball B of M" can also be 
defined by applying the Lorentz boost 



B{u) 



7u 7uU" 
7uu I + j^uu^ 



/ 1 \ 

to and obtaining 



V 



Bin) 

where t = 7u(l + u'^v). 

If two vectors u and v are parallel, that is, there exists a nonzero constant A such that 
V = Au, then 

u -|- V 

ue V 



1 + u^v ■ 

To capture abstractly Einstein addition of velocities in special relativity, A. A. Ungar has 
introduced and studied in several papers and books structures that he has called gyrogroups; 
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see [5] and its bibliography. His axioms are reminiscent of those for a group, but gyrogroup 
operations are nonassociative in general. 

Definition 2.3. A triple {G, ©, 0) is a gyrogroup if the following axioms are satisfied for all 
a,b,c £ G. 

(Gl) 0©a = a©0 = a (existence of identity); 

(G2) a © {—a) = (—a) © a = (existence of inverses); 

(G3) There is an automorphism gyr[a, 6] : G — )• G for each a,b €z G such that 



A gyrogroup is uniquely 2-divisible if for every b G G, there exists a unique a G G such that 



We call gyr[a, 6] the gyroautomorphism or Thomas gyration generated by a and b. 

Remark 2.4. It has been shown that gyrocommutative gyrogroups are equivalent to Bruck 
loops jl], i.e, a loop is a Bruck loop if and only if it is a gyrocommutative gyrogroup 
with respect to the same operation. It follows that uniquely 2-divisible gyrocommutative 
gyrogroups are equivalent to S-loops, uniquely 2-divisible Bruck loops. 

A. A. Ungar has shown in [5, Chapter 3] by computer algebra that Einstein addition 
on the ball B is a gyrocommutative gyrogroup operation, and the gyroautomorphisms are 
orthogonal transformations preserving the Euclidean inner product and the inherited norm. 
We call (B,©) the Einstein (gyrocommutative) gyrogroup, where © is defined by ()2.ip . 

In references [5], [6], and [7] the Einstein vector addition is mostly defined on the open 

unit ball B of M". It turns our interest to extending the Einstein addition on the closed 

unit ball B. Substituting Ou = — in the definition (12. Ih we have an alternative expression 

7u 

of the Einstein vector addition 



This is well defined for all (u, v) G B x B excluding the antipodal vectors. Since u© v = u 
for any v G B if ||u|| = 1, we are able to define the Einstein addition naturally for the 



a © (6 © c) = (a © 6) © gyr[a, b]c (gyroassociativity) ; 

(G4) gyr[0,a] = id; 

(G5) gyr[a © 6, 6] = gyr[a, b] (loop property). 
A gyrogroup {G, ©) is gyrocommutative if it satisfies 

a © 6 = gyr[a, b]{b © a) (gyrocommutativity). 



a © a = 6. 
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antipodal vectors by continuity such that 

U © V = u 

for any v € B whenever ||u|| = 1. Hence, the extended Einstein addition of u and v in B 
can be defined as 

' (121]) if u E B 
u if Hull = 1 



u 



Remark 2.5. The closed unit ball B with respect to the extended Einstein addition is a 
binary system, but not a gyrogroup since it has no unique inverse. 

3. Block vectors and density matrices 

A qubit is a member of a 2-dimensional Hilbert space, or a two-state quantum system. 
A qubit density matrix is a 2-by-2 positive semidefinite Hermitian matrix with trace 1. 
Indeed, any 2-by-2 Hermitian matrix of trace 1 must have a parametrization of the form 

_ 1 / 1 + W3 V\-iV2\ 

'^\v\^iv2 1 - f 3 / ' 

where 

V = V2 G M^. 
Us / 

So the qubit density matrix can be described as pv for some v E such that ||v|| < 1. In 
this case the vector v is known as the Bloch vector or Bloch vector representation of the 
state /?v 

Remark 3.1. Via the characteristic polynomial of the qubit pvj we obtain that its eigen- 
values are 

1-F||v|| 1 - ||v|| 



|2 



2 ' 2 

and its determinant is 

l-||v„ 
det yOv = ^ • 

So the mixed states are parameterized by the open unit ball B in M"^. 
The Pauli matrices are given by 
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where i = \/— 1. The parameterization of qubit density matrices can be written in an 
alternative method using the Pauh matrices: 



where y'^ a = viax + V2(Jy + V2,(Tz for 





( V, \ 




/ 




\ 


V = 








ay 














I 



Let ID) = {pu '■ u S B} be the set of all qubit density matrices representing the mixed 
states. We define a binary map : ID) x D — )- B as 

1 



PVL /Ov 



tr (pu * Pv 



-Pu * Pv 



1/2 1/2 

where /?u * Pv = Pu PvPu • 



Remark 3.2. From the fact that \x AB = \xBA for any matrices A and -B, we have [U 
Eq. (9.67)] 



tr(pu*Pv) = trpupv 



1 + u^v 



This gives us that the binary map is well defined since 1 + u^v 7^ whenever u, v E B. 

Since every element in D is a Hermitian positive definite matrix, it has a unique square 
root. We have the explicit form and can prove it via the direct matrix computation of 
squaring the form (j3.3p . 



Lemma 3.3. For any v S B the square root of the qubit density matrix is uniquely given 
by 



1/2 

Pv 



7v 



l + 7v 



(3.3) 



We now see an isomorphism between the open unit ball B with the Einstein velocity 
addition © and the binary system (D, 0). 

Theorem 3.4. The map p : (B,0) — )• (B,0), v 1— t- pv is an isomorphism. 

Proof. Obviously the map p is a bijection. We need to show that Pu®v = Pu Pv for any 



u 



Ul U2 U3 ,V 



Vl V2 V3 



G B. 
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We set X := . By Lemma 13.31 we have 

1 + 7u 



1/2 _ 1 ^ ( ^3 + ^ Ul- 



^ \ Ul + IU2 -Us + - J 

It is enough to check the (1,1) and (1,2) entries of pu Pv since pu Pv is Hermitian. 
Let us first compute the (1,1) entry of pu * /^v Then 

^3 + - ) (1 + V3) + (^^1 - iu2){vi + iV2) I ( tis + - 



f { (^^ ^ x) ~ ^^^^ ^ '■^^ ~ iu2){l - f3)| {ui + iU2) 
us + (1 + vs) + 2 (^U3 + (nivi + U2V2) + {uj + ui)(l - vs) 

"|u||2 + 2(u^v)n3 - \\ufvs) + -us + -(u^v) + ^(1 + Vs) 
|||u||2(l - vs) + ^(u^v)n3 + ^us + ^(u^v) + ^(1 + vs) 

^ ^"^3 +777— V U3 + -U3 + -(U V + — I+V3) 

87u 4(1 + 7u) 4 4 87u 

^ J (1 + u^v) + n3 + — ^3 + -^^(u^v)m3 I . 

7u 1 + 7u J 

By Remark 13.21 the (1,1) entry of p^ Pv is 

7 I (1 + U^v) +US + —Vs + — ^(u'^v)ti3 



1 + vy 4 [ 7u 1 + 7u 

' ' + 1 ^ T "3 + — ^^3 + -— (U 

1 + u-* V V 7u 1 + 7u 



= ^(1 + (U©V)3), 

where (u © v)j represents the j th coordinate of Einstein vector addition u © v. 
Similarly, let us compute the (1,2) entry of pu * Pv Then we have 

I I ^^i3 + (1 + Vs) + {ui - iu2){vi + ro2)| {ui - iU2) 

+ I I flt3 + ^ j (t'l - 'i'V2) + [ui - iU2){l - Vs) \ i-Us + ^ 
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The real part of the (1,2) entry is 

I I ^^^3 + (1 + V3)ui + {uj - ul)vi + 2uiU2V2 

+ ^ I f ^^3 + H ^^1 (-U3 + ^) + ~ + ^ 



x[2 / T \ iiii2 ^ 

o \ -ui + 2(u v)ui - u ui + —vi 
8 [x x^ 



and the imaginary part is 



X ( ( 1 \ 

Us + - (1 + t'3)n2 + (n? - ti|)f2 - 2u\U2V\ 

X J 

^ <j -2n2 ( n3^;3 + i j + (uf - n| + n^)t;2 - 2uiU2fi - -^^2 



x[2 ( T \ iiii2 

- < U2 - 2(U V)n2 + U V2 

8 X x"^ 



= -7 ""2 H f2 + — (U V)U2 

By Remark 13.21 the (1,2) entry of pu Pv is 

2 \ 1 + u^v,, ... . . 1 



1 + u^v/ 4 ■((u©v)i -i(u© v)2) = -((u©v)i -i(u© v)2). 
Therefore, we conclude that O Eq. (9.23)], 

\( l + (u©v)3 (u© v)i -i(u© v)2 \ 

2 y (u© v)i +z(u© v)2 l-(u©v)3 j 

for any u, v € B. This means the map /3 is a homomorphism. □ 

Remark 3.5. From Theorem 13.41 we obtain that D is a gyrocommutative gyrogroup with 
respect to the operation defined by 

^ 1 1/2 1/2 

Pu Pv = — -T-TFi 1/2^^" P^P^ 

tr pu PvPu 
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On B, moreover, the identity is (1/2)/ and the inverse for is 

P-u = ^Pu'- (3.4) 

47u 

4. The rapidity metric 

The Einstein gyrogroup on the open unit ball B admits an analytic scalar multiplication 
given by 

(t,u) H> t • u = tanh(ttanh"^ ||u||)(u/||u||) (4.5) 

for u 7^ 0, and for t = or u = (see [S] Chapter 6]). Expressing the magnitude of the 
velocity parameter u by the hyperbolic parameter (/)u called the rapidity of u G B, 

(pu = tanh""*" ||u||, 

we have t ■ u = tanh(t(/)u)(u/||u||), or (pt-u = t4>M- 

For the Einstein gyrogroup (B,©,0), A. A. Ungar considers what we call the Ungar 
gyrometric defined by q(u, v) = || — u© v||. He also defines what we call the rapidity metric 
by d(u, v) = tanh^^ f)(u,v). It is known as the Cayley-Klein metric on the Beltrami-Klein 
model of hyperbolic geometry (see [2i] ) , or the Bergman metric on the symmetric structure 
B with symmetries S'w(v) = u© (— v) for some u = 2 • w (see |8j). For real numbers s,t, 
we define 

s + t 

s © t 



l + st' 

the restricted Einstein addition analogous to the Einstein sum of parallel vectors. We see 
some properties of rapidity metric on B. 

Lemma 4.1. The following properties hold for all u, v, w G B. 

(i) 0< £)(u,v),d(u,v) 

(ii) q{u, v) = <^ d{u, v) = 44> u = V 

(iii) £>(u,v) =g{^r,u), d{u,^r) = d(v, u) 

(iv) ||u © v|| < ||u|| © ||v|| q{u, w) < g{u, v) © q{v, w) d{u, w) < d(u, v) + d{v, w) 

(v) q{u © V, u © w) = Q{\r, w) and d(u © v, u © w) = d(v, w) 

(vi) d{0,r- w) = |r|d(0, w). 

We establish the relation of metrics on gyrocommutative gyrogroups under an injective 
homomorphism. 



Lemma 4.2. Let (Gi,©,0) and (G2,©,0) be gyrocommutative gyrogroups each equipped 
with a metric invariant under left translations. If ip : Gi ^ G2 is an injective gyrogroup 
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homomorphism and if there exists k > such that (i(0,^(u)) > Kd(0,u) for each u G Gi, 
then d{ip{u),'ip{\-)) > Kd{u,v) for all u, v € Gi. 

Proof. Let u, v € Gi and set w = — u© v € Gi. By hypothesis d{0,^{w)) > Kd{0,w). By 
invariance of the metrics under left translation d{0, w) = d{u © 0, u © (— u © v)) = fi(u, v) 
and similarly 

d(V'(0), V(w)) = (i(^(u), V(u) © V'l-u © v)) = d(V'(u), V(v)). 
We conclude that d(^(u), V'(v)) = (i(0,V'(w)) > Kd{0,w) = Kd(u,v). □ 

5. The trace metric for qubit density matrices 

Here we review the definition of the trace metric on the open convex cone of (Hermitian) 
positive definite matrices. The trace metric on Q is determined locally at the point A by 
the differential 

ds = \\A-'/^dAA-'/^F, 

where || • \\f means the Frobenius or Hilbert-Schmidt norm. This relation is a mnemonic 
for computing the length of a differentiable path 7 : [a, 6] — )• 

L(7)= [''h-'/\th'{t)r'^\t)\\pdt. 

J a 

Based on the above notion of length, we define the trace metric 6 between two points A 
and S in 17 as the infimum of lengths of curves connecting them. That is, 

6{A,B) := inf{ ^(7) : 7 is a differentiable path from A to B }. 

R. Bhatia and J. Holbrook have established some basic properties of the trace metric on 
n in [1]. 

Lemma 5.1. For any A,B & Q and any invertible matrix X, 

(i) 5{T x{-^)^^ x{B)) = 6{A,B), where Tx is a congruence transformation by X. 

(ii) 5{A,B) = \\\ogA-^''^BA-'^/'^\\F = \\\og{A-^B)\\F. 

(iii) 5{A^'\B^'^)<]^6{A,B). 

We now consider the trace metric on the set D of all invertible qubit density matrices, 
and then see the connection with the rapidity metric on the open unit ball B. 

Proposition 5.2. For any u, v € B, 

«(.„.P.)<{ln^O+..^(f)}'^ 
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where x = , a = -—^ , and b = -— r • 

2 472' 472 

Proof. By the equation ()3.4p and Theorem 13.41 we have 

-1/2 -1/2 x2 1/2 1/2 

Pu PvPu = (27u) P-uPvP-u 

2 /I - U^v\ 

= 47u ( 2 j ^-"®^- 

We set K = 47^ ^ ^ ^ andp= -(1 + || —u©v||), wherep and 1—p are the eigenvalues 

of/9_uev Then 

^2 111 -1/2 -1/2||2 

= 2 In^ K + 2 In K{lnp + In (1 - p)} + In^ p + hi^ (1 - p) 
< 2 In^ K + 2 In X In - p)) + In^ - p)) 
= ln'^K + {In K + In - p))}^ 
= In^i^ + ln2 {Kp{l-p)) 



The inequahty holds since Inp < whenever < p < 1. The last equality follows from the 
property 7u®v = 7u7v(l + u^v). Indeed, 

Kp{l-p) = A^i 



2 I \ 4 
l-u'^v\ 1 



2 J 47I,, 
1 



2 

V 



272(1 - U^V) 
X 

b' 



□ 



Remark 5.3. From Remark 13.21 and Remark 13. H we see that 

1 - u^v 

X = 2 = (/^-uPv) 

a = -—5- = det Pu and b = — ^ = det pv 

472 474 
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Theorem 5.4. For any u, v G B, 

d(u,v) < --^(5(pu,Pv)- 

Proof. We have seen in Remark 13.51 that (D, 0) is a gyrocommutative gyrogroup via an 
isomorphism 

p : (B,e) ^ (B,0). 

By Lemma [4. If v) and Lemma [5. If i) we have the rapidity metric on B and the trace metric 
on B are both invariant under left translations. By Lemma |5. If ii) and Remark 13.11 

6{{l/2)I,p^f = \\log{2p^)\\l 

= ln2(l + ||u||)+ln2(l-||u||) 

>lln^l + HI 



2 l-||u|| 
= ^(2tanh-i ||u||)2 = 2d{0,uf. 

So (5((l/2)/, pu) > \/2fi(0,u). Therefore, by Lemma l42] we conclude 

S{pu,Pv)> \/2d(u,v). 
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